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SUMMARY 

Data  from  eighteen  reported  Wind  Tunnel  and  Free  Flight  tests 
on  heat  transfer  to  blunt  cones  and  hemispheres  are  correlated  with  simple 
theoretical  formulae.  The  Mach  number  and  Reynolds  number  ranges  for  these 
tests  are  from  1.33  to  10.2  and  8  x  10^  "to  4*5  x  10^  respectively.  The 
blunt  cone  models  have  semi-angles  between  4*75°  and  15°  and  nose  to  base 
radius  ratios  between  0.125  and  0»8. 

For  the  Free  Flight  tests,  estimates  of  the  error  in  measured 
heat  transfer  parameters  are  made  and  only  points  showing  errors  of  less 
than  1 O^b  are  retained  in  the  analysis.  The  Wind  Tunnel  results  are  claimed 
to  have  errors  within  this  limit. 

In  this  report  the  Nusselt  number  -  Reynolds  number  relationships 
are  investigated  and  Reynolds  Analogy  Factors  determined  for  a  selection 
of  the  tests.  In  these  correlations  corrections  for  compressibility  effects 
are  made.  The  ratios  of  measured  Nusselt  numbers  to  the  corresponding 
theoretical  imcompressible  Nusselt  numbers  are  also  determined.  The 
principal  conclusions  made  are  as  follows 

1 )  For  laminar  and  turbulent  flow  over  hemispheres  and 
hemispherical  nose  caps,  the  Nusselt  and  Reynolds  numbers  are 
related  by  simple  power  law  formulae. 

2)  For  laminar  flow  Wind  Tunnel  tests  on  conical  skirts,  the 
Nusselt  numbers  lie  between  flat  plate  and  sharp  cone  theories 
when  plotted  against  Reynolds  number.  As  the  value  of  the  ratio 
of  nose  radius  to  base  radius  increases,  Nusselt  numbers  move 
away  from  the  values  predicted  by  sharp  cone  theory. 

3)  For  Free  Flight  tests  on  conical  skirts  with  laminar  flow,  the 
Nusselt  and  Reynolds  numbers  are  best  correlated  by  sharp  cone  theory 
when  they  are  based  on  the  length  measured  from  the  virtual  apex 

of  the  cone  skirt. 
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For  both  Wind  Tunnel  and  Free  Flight  tests  on  conical  skirts 
with  turbulent  flow,  the  Nusselt  and  Reynolds  numbers,  based  on  the 
lengths  measured  round  the  model  surfaces,  are  related  by  the 
four-fifths  power  law  for  flat  plates. 

Correlations  of  Free  Flight  measurements  on  conical  skirts 
are  independent  of  whether  the  local  flow  conditions  are  based  on 
expansion  from  stagnation  pressures  corresponding  to  normal,  or  to 
conical  shock  wave  conditions. 

The  ratio  of  experimental  Nusselt  numbers  to  theoretical 
incompressible  Nusselt  numbers,  when  multiplied  by  a  factor  accounting 
for  compressibility  effects,  are  independent  of  the  local  values 
of  Mach  number,  Reynolds  number  and  Wall  temperature  ratio. 

Reynolds  Analogy  Factors  appear  to  be  independent  of  local 
Mach  number  and  Reynolds  number. 

Wall  temperatures,  on  thin  skin  Free  Flight  models,  can 
generally  be  predicted  to  an  accuracy  of  10^  using  expressions  for 
heat  transfer  coefficients  given  by  Monaghan  and  Hill's  wall 
temperature  equations. 
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INTRODUCTION 

Hypersonic  flight  is  inevitably  associated,  with  large  heat  transfer 
rates,  particularly  round  the  forward  regions  of  Missiles  and  Space¬ 
craft.  This  problem  has  been  partly  solved  by  the  introduction  of 
blunted  noses  on  the  vehicles  and  by  the  introduction,  for  re-entry 
vehicles,  of  ablating  heat  shields.  The  problem  of  predicting  the 
convective  heat  transfer  distribution  around  blunt  cones  and  hemispheres 
has  therefore  been  extensively  investigated,  both  theoretically  and 
experimentally,  in  recent  years. 

Broadly,  one  may  divide  these  investigations  into  two  classes; 

(a)  heat  transfer  at  Mach  numbers  greater  than  ten  and  (b)  heat 
transfer  at  Mach  numbers  less  than  ten.  In  the  first  category,  real 
gas  effects,  such  as  dissociation,  dominate  the  problem  and  the 
'cool  wall'  approximation  is  valid.  In  the  second,  real  gas  effects 
are  minor  but  the  specific  heats  are  variable.  The  'cool  wall' 
approximation  is  not  valid  in  this  second  category  of  problems. 

This  report  presents  the  results  of  an  investigation  into  the 
latter  class  of  heat  transfer  problems  just  mentioned.  The 
investigation  has  been  restricted  to  hemispheres  and  homispherically 
blunted  conos,  which  probably  represent  the  most  extensively  tested 
bluff  body  shapes.  Although  no  claim  is  made  that  all  the  experimental 
evidence  has  been  investigated,  it  is  felt  that  a  sufficiently  large 
number  of  reported  tests  have  been  analysed  for  the  conclusions 
drawn rto  be  representative. 

For  both  Wind  Tunnel  and  Free  Flight  tests  the  data  from  eighteen 
reported  experiments  on  hemispheres  and  blunted  cones  has  been 
analysed  with  the  object  of  determining,  for  laminar  and  turbulent 
flow,  the  followings- 

1 )  The  variation  of  Nusselt  number  with  Reynolds  number  and  the 
optimum  lengths  on  which  to  base  these  parameters. 

2)  The  effects  of  local  Mach  number,  wall  temperature  ratio  and 
Reynolds  number  on  the  ratios  of  the  compressible  Nusselt  numbers 
obtained  by  experiment,  to  the  incompressible  theoretical  estimates. 

3)  For  conical  skirts;  the  effect  of  basing  the  local  flow  conditions 
on  isentropic  expansion  from  either  normal  or  conical  shock  wave 
conditions. 

4)  The  effect  of  nose  to  base  radius  ratio  on  the  Nusselt  number  for 
conical  skirts. 

5)  Evidence  of  Reynolds  Analogy  for  a  selection  of  the  experiments 
analysed. 
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Comparison  of  theoretical  and  measured  wall  temperatures 
on  Free  Flight  Models  with  thin  skins. 


Section  3  gives  details  of  the  methods  used  to  reduce  and  present  the 
data,  and  in  Section  4  correlated  results  are  discussed  in  detail.  The 
final  conclusions  drawn  from  the  investigation  are  presented  in  Section  5* 
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2.  Notation 
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c 

Cf 

CP 

CP. 

G 

h 

h' 

i 

^-r 

k 

M 

Nu 

Nue 

%a 

Nut 


Nu 


LD 


roN 

Pr 


Specific  heat  of  skin  material 

Local  skin  friction  coefficiont 

Specific  heat  of  air  at  constant  pressure 

Pressure  coefficient 

Thermal  capacity  of  skin  p  a;'x  * 

r  8  8  w‘ 

Heat  transfer  coefficient  based  on 

temperaturo 

Heat  transfer  coefficiont  based  on  enthalpy 
Enthalpy 

Recovery  enthalpy 
Thermal  conductivity  of  air 

Mach  number 
Nusselt  number 

Nusselt  number  based  on  s 

Nusselt  number  based  on  sa 
Nusselt  number  based  on  s^ 

Nusselt  number  based  on  D 
Static  pressure 

Total  pressure  behind  a  normal  shock  wave 
Prandtl  number 

Convective  heat  transfer  rate 


CHU.lb  1  °K  1 


CHU.lb  1  °K_1 


CHU.ft"2  °K  1 

CHU.ft"2  sec-1  °K"1 
lb.  ft-2  sec-1 


-1 


CHU.  lb 


1  —1  n  —1 

CHU.ft"  sec.  K 


hs/Kl  or  h'  Cp^/l^ 
h'  Cpi  s/K1 
h»  Cpi  Sa/K1 
h'  Cpi  St/K1 
hD/K1 


lb.  ft 


-2  / 


-2 


lb.  ft 


CHU. ft-2. sec"1 


Recovery  factors  Local  body  radius  (Fig. l) 

Nose  radius  for  blunted  cone  or  hemisphere  (Fig. 1.)  ft. 

Base  radius  of  blunted  cone  (Fig. 1.)  ft. 

,  Reynolds  number 

js  Reynolds  number  based  on  s  p  Uis/ 

Ren  Reynolds  number  based  on  s-  o  tt  / 

a  a  Pi  Ui3a/ii1 

Ret  Reynolds  number  based  on  p  U  s  ^ 

s  Surface  distance  measured  from  stagnation  point  (Fig. 1)  ft. 
Sr,  Distance  measured  from  virtual  apex  of  blunt  cone  (Fig.l)  ft. 


r 

R 

Rb 

Re 

R 


st 

S 

Sip 

t 


Distance  measured  from  transition  point  (Fig.  1 ) 

Reynolds  Analogy  Factor 

Stanton  number 

Time 


ft. 


h/p  U  C  or  h' 

/ri  1  p 


P 
sec. 
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T 

Tr 

To 

U 

P 

r 

8 


Temperature 

Recovery  temperature 

Stagnation  temperature 

Velocity  along  s  co-ordinate 

Stagnation  point  velocity  gradient 

Ratio  of  specific  heats  for  air 

Time  interval  in  wall  temperature 

calculation 


°K  (except  where  otherwise  stated) 

°K 

°K 

(ft-sec^ ) 

-1 

sec 

sec. 


e 

9c 

P 

v 

P 

c r 


t 


Surface  emissivity 
Cone  semi-angle 
Dynamic  viscosity 
Kinematic  viscosity 
Density 


degrees 

lb.  ft-1  sec-1 

ft.2sec.  1 


lb.  ft 


-3 


Deviation  in  a  parameters  also  ^  q  . 

Stefan-Boltzmann  radiation  constant  CHU.ft  sec.  K  ^ 


Skin  thickness 

Angular  co-ordinate 

Pressure  gradient  parameter 
defined  in  Appendix  VI 


ft. 

degrees 


Subscripts 

o  Stagnation  conditions 

1  Edge  of  boundary  layer 

oo  Free  stream  conditions 

c  Compressible  flow 

i  Incompressible  flow 

f.p.  Plat  plate  conditions 

w  Wall  conditions 

Superscript 

*  Properties  based  on  intermediate  enthalpy  conditions 


N.B.  All  other  symbols  used  in  the  report  are  defined 
in  the  text. 
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3.  Data  Reduction  and  Presentation  of  Results 

3. 1 .  Data  Reduction  Programme 


Experimental  data  on  heat  transfer  to  blunted  cones  and 
hemispheres  from  eighteen  sources  (Refs.  1-18)  has  been  analysed, 
and,  since  the  results  in  these  reports  are  presented  in  terms  of 
various  parameters,  a  dr.ta,  reduction  programme  was  prepared  to 
reduce  all  the  data  to  common  sets  of  non-dimensional  parameters. 
Wherever  possible  the  local  flow  conditions  were  calculated 
from  the  isentropic  flow  relations  using  experimental  pressure 
distributions.  In  cases  where  the  latter  were  unavailable,  modified 
Newtonian  theory  was  used.  Referring  to  Fig.  1 ,  the  Newtonian 
pressure  coefficient  distribution  is  given  by 


C 


P 


C  .  cos2  y 

y  max 


.0) 


where  CL  max  =  Pnw  -  p  <» 

p  _2£! - -  - (2) 

i  .yPoo  Mm  2 

and  P  „  was  estimated  from  the  normal  shock 
oN 

wave  relations  withy  ■=  1.4* 


On  the  conical  skirts  of  slightly  blunted  cones,  where  the 
pressures  were  not  measured,  the  pressure  was  assumed  to  take  the 
conical  value  given  by  Kopal  (Ref. 19)» 

The  Eckert  intermediate  enthalpy  technique  (Ref.  20 )  was  used 
throughout  the  correlations,  the  relevant  properties  being  evaluated 
at  a  temperature,  T,*  derived  from  the  enthalpy,  i,*  which  is 
given  by 

i*  =  i1  +  0.5  (iw  -  i.,)  +  0.22  (ir  -  i.,) 

_ (3) 

with  ir  =  i1  (1  +  r  ( y  -  1 )  M  2  ) 

2 


where  r  =  O.85  for  laminar  flow 
and  r  =  O.89  for  turbulent  flow. 

#  1 

(NOTE:-  For  laminar  flow,  strictly  r  *=  Pr  2  but  r  «=  O.85  is 

well  within  the  experimental  accuracy  and  simplifies  the 
computing) 
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As  in  the  analysis  carried  out  in  Ref.  21 ,  only  experimental 
results  having  errors  of  10%  or  less  were  retained  in  the  present 
analysis.  For  tho  Free  Flight  tests,  the  errors  in  measured 
Stanton  numbers  wore  estimated  on  the  same  basis  as  that  used  in 
Refs.  21  and  22.  This  implies  that  the  errors  in  Stanton  number 
due  to  errors  in  tho  estimated  pressure  distributions  havo  been 
ignored.  Tho  equation  for  estimating  the  errors  together  with 
detailed  sources  of  error  for  Free  Flight  tests,  given  in  Ref.  22, 
are  summarised  in  Appendix  I.  Errors  arising  from  conduction 
along  tho  skin  close  to  the  stagnation  point  were  small  inrcases 
where  it  was  possible  to  assess  these.  A  typical  case  is  discussed 
in  Section  4. 3.1., the  relevant  equations  being  derived  in  Appendix  II. 

3.2.  Presentation  of  Results 

Details  of  all  the  tests  analysed  are  presented  in  Fig.  2,  which 
also  includes  a  key  to  the  symbols  used  in  the  graphs.  As  far  as 
possible  each  source  of  data  has  been  given  the  same  symbol  through¬ 
out  the  correlations. 

It  has  been  found  convenient  to  present  most  of  the  results  in 
groups  as  follows :- 

(i)  Results  for  hemispheres  and  nose  caps  of  blunted  cones 
with  laminar  flow. 

(ii)  As  for  (i)  with  turbulent  flow. 

(iii)  Results  for  cone-skirts  with  laminar  flow. 

(iv)  As  for  (iii)  with  turbulent  flow. 

In  addition,  where  sufficiently  large. .numbers  of  both  Froo  Flight 
and  Wind  Tunnel  results  are  available,  these  are  presented 
separately  for  clarity. 

3.2.I.  Variation  of  Nusselt  number  with  Reynolds  number 

*  4- 

The  parameter  Nu  Pr  3  ^determined  for  hemispheres, 
hemispherical  caps  and  cone  skirts,  is  plotted  in  Figs  3-13 
against  Re*  based  on  various  lengths.  (See  Fig.  1). 

Figs.  3  and  4  present  Wind  Tunnel  and  Free  Flight 
results  respectively,  for  hemispheres  and  hemispherical 
caps  with  laminar  flow,  tho  abscissa  in  both  cases  being 
Re  based  on  s.  In  Figs.  5  and  6,  Wind  Tunnel  and  Free 

Flight  results  for  hemispheres  and  hemispherical  caps 

* 

with  turbulent  flow  are  plotted  against  Re  based  on 
st.  (Fig.  6.) 
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*  *  —  1 

Values  of  Nu  Pr  ^  for  conical  skirts  with  laminar  flow  are 
presented  in  Figs.  7-10;  Figs  7  and  9  show  plots  of  Wind  Tunnel  and 
Free  Flight  results  respectively,  against  Re*  based  on  s,  whilst 
Figs.  8  and  10  show  these  same  results  plotted  against  Re*  based 
on  s„  •  Figs  11-13,  which  are  the  final  graphs  in  this  first  group, 
present  Free  Flight  and  Wind  Tunnel  results  for  conical  skirts  with 
turbulent  flow,  tho  abscissae  being  Re  based  on  s,  sa  and  s^. 

respectively. 

In  Figs.  9-13,  the  full  symbols  correspond  to  local  flow 
conditions  based  on  flow  through  a  conical  shock  wave  originating 
from  the  vertex  of  a  cone  having  the  same  apex  angle  as  the  cone 

skirt  under  consideration.  The  remaining  symbols  correspond  to 
local  flow  conditions  based  on  isentropic  flow  from  stagnation 
conditions  behind  the  normal  section  of  the  bow  shock  wave.  This 
is  essentially  the  method  suggested  by  Moeckel  in  Ref.  23. 

The  equations  of  the  theoretical  curves  shown  in  Figs.  3-13 
are  derived  in  Appendix  III,  and  on  each  figure,  deviations  of 
+  50$  from  the  appropriate  theoretical  curve  are  shown  in  order 
to  indicate  the  spread  of  the  data. 

3.2.2.  Variation  of  the  Laminar  Heat  Transfer  Parameter 
WuD/(  p  D2/  v  with  s/R 

Wind  tunnel  values  of  the  laminar  heat  transfer 
2  i 

parameter,  Nup/(pD  /v  )  for  hemispheres  and  hemispherical 
caps  are  presented  in  Fig.  14  as  a  plot  against  s/R^the 
normalised  distance  along  the  body  surface.  Since  it  depends 
upon  the  local  Mach  number  and  the  local  wall  temperature 
ratio,  theoretically  this  parameter  does  not  collapse  the 
results  on  to  a  single  cuive.  However,  this  parameter  is 
often  used  to  present  Wind  Tunnel  results  of  laminar  heat- 
transfer  measurements,  and  has  been  included  for 
completeness  in  the  present  report.  For  the  sake  of  comparison 
the  theoretical  curve,  obtained  by  Korobkin  in  Ref.  31, 
is  shown  in  Fig.  14*  It  should  be  noted  here  that  the 
corresponding  turbulent  heat  transfer  parameter  also  fails 
to  collapse  the  results  on  to  a  single  curve  when  plotted 
against  s/R. 

3.2.3.  Variation  of  tho  Ratio  of  Compressible  to  Incompressible 
Nusselt  numbers  with  s/R 

As  stated  above,  the  parameter  Nujj  J  (pD2/v)2  does  not 
collapse  the  results  when  plotted  against  s/R.  However,  the 
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parameter  (Nug/Nu^)  (T*/T-|)a  should  (as  shown  in  Appendix  IV) 
collapse  the  data  for  hemispheres  on  to  a  straight  line  when 
plotted  against  s/R,  with  a  =  0.13  for  laminar  flow  and  0.64 
for  turbulent  flow.  By  considering  this  particular  parameter, 
the  effects  of  local  Mach  number  and  wall  temperature  ratio 
should  be  eliminated. 

In  Figs.  15  and  16,  Wind  Tunnel  and  Free  Flight  results 
respectively,  are  presented  for  hemispheres  and  hemispherical 
eaps  with  laminar  flow,  and  in  Fig.  18  results  obtained  from 
Wind  Tunnel  and  Free  Flight  tests  for  the  above  bodies  with 
turbulent  flow  are  presented.  In  Figs.  15>  16  and  18  deviation 
of  -  50fc  from  the  theoretical  line  are  again  shown  in  order  to 
indicate  the  spread  of  the  data. 

Wind  Tunnel  values  for  conical  skirts  with  laminar  flow  are 
shown  in  Fig.  1%  whilst  laminar  results  on  conical  skirts 
obtained  from  Free  Flight  tests  are  shown  in  Figs.  20  and  21. 

In  Fig.  20,  all  the  local  flow  conditions  are  based  on  the 
assumption  that  the  surface  streamlines  have  passed  through  a 
normal  shock  wave,  whereas  in  Fig.  21  it  has  been  assumed  that 
the  surface  streamlines  have  passed  through  a  conical  shock 
wave  (c.f.  Section  3.2. 1.)  Fig.  22  shows  a  plot  of 
(NusAui)(T*/Ti)°‘64  against  s/R  for  conical  skirts  with 
turbulent  flow,  the  results  being  obtained  from  Wind  Tunnel 
and  Free  Flight  tests.  Here,  both  normal  and  conical  shock 
wave  conditions  have  been  used  to  obtain  two  sets  of  values 
of  the  above  parameter;  both  are  plotted.  For  cone  skirts  this 
parameter  does  not  collapse  the  results  on  to  a  single  line, 
but  as  explained  in  Appendix  IV,  the  results  should  be 
between  two  lines  representing  the  flat  plate  and  conical 
limits  of  the  parameter. 

In  Fig.  17,  Nus/Re  ^  is  plotted  against  s/R  for  the 
s 

results  of  Reference  6,  with  and  without  heat  conduction 
effects  accounted  for  in  the  region  near  the  stagnation  point. 
The  derivation  of  the  equation  for  estimating  the  effect  of 
longitudinal  heat  conduction  is  presented  in  Appendix  II. 
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3.2.4.  Reynolds  Analogy  Factor 

Reynolds  Analogy  Factor  was  estimated  in  two  ways,  both 

involving  the  use  of  experimental  Stanton  number 

measurements.  (See  Appendix  V)  The  first  was  the  evaluation 

of  the  skin  friction  coefficient  from  a  flat  plate  formula, 

and  the  second  was  the  correction  of  this  coefficient  for 

pressure  gradient  effects  as  described  in  Appendix  VI. 

The  function  ljr ,  derived  in  Appendix  VI,  which  describes 

the  pressure  gradient  effects,  is  shown  in  Fig.  23  for  a 

hemisphere-cone  of  4*75°  semi-angle  with  laminar  flow  at 

a  free  stream  Mach  number  of  3.12  and  Tw  ratio  of  1.0. 

To 

Results  for  the  two  methods  of  estimating  Reynolds 
Analogy  Factor  for  a  selected  number  of  Wind  Tunnel  tests 
on  hemispheres  and  hemispherical  caps  are  presented  in 
Figs.  24-27.  Figs  24  and  26,  which  show  plots  of  S  against 
local  Mach  number  and  Reynolds  number  respectively, 
also  include  the  boundaries  within  which  lie  the  points 
corrected  for  pressure  gradient  effects.  Values  of 
Reynolds  Analogy  Factor  with  the  correction  for  pressure 
gradient  effects  are  presented  in  Figs.  25  and  27  as 
plots  against  local  Mach  number  and  Reynolds  number 
respectively. 

In  Figs.  28-33*  the  estimates  of  Reynolds  Analogy 
Factor  on  the  conical  skirts  are  presented,  the  values 
being  obtained  from  the  equations  derived  in  Appendix  V. 
These  estimates  are  plotted  against  both  local  Mach  number 
and  local  Reynolds  number.  Again,  the  full  symbols  of 
Figs.  30,31,32  and  33  correspond  to  local  flow  conditions 

based  on  flow  through  a  conical  shock  wave.  On  all  these 

_  2 

plots  the  value  S  -  Pr  3  ,  suggested  by  Colburn  (Ref. 24), 
is  shown  for  the  purpose  of  comparison. 

3.2.5*  Wall  temperature  measurements 

Wall  temperature  results  for  four  selected  Free  Flight 
tests  are  shown  in  Fig.  34.  The  results  are  presented  as 
the  ratio  of  measured  to  estimated  wall  temperature 
plotted  against  a  reduced  time  of  flight,  t,  defined  as 

t  «=  (t  -  ■fcmin)/(l;max  -  't'min) 

where  t  min  =  initial  time  for  data  recording 
and  t  max  =  time  at  which  maximum  measured  wall 
temperature  occurs. 

The  method  used  for  estimating  the  wall  temperatures 
is  described  in  Appendix  VII. 
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4.1.  Variation  of  Nusselt  number  with  Reynolds  number 

4.1.1.  Results  for  hemispheres  and  nose  caps  (Figs  3-6) 

Both  Wind  Tunnel  and  Free  Flight  data  for  hemispheres 
with  laminar  flow  correlate  reasonably  well  with  the  simple 
hemisphere  theory,  as  seen  in  Figs  3  and  4*  The  low 
points  from  Ref.  8,  at  the  high  Reynolds  number  end  of 
Fig.  3,  coming  from  one  source,  cannot  be  regarded  as 
strong  evidence  for  modifying  the  theory  in  this  region. 

It  would  Appear  that  there  is  some  consistent  error  in 
the  measurements,  possibly  due  to  conduction  effects. 

(As  shown  in  Appendix  II,  conduction  may  lead  to  an 
under-estimate  of  the  Nusselt  number).  If  these  points 
are  discarded  then  the  theory  of  Ref.  25,  is  well 
verified  by  the  experimental  results  in  the  Reynolds 
number  range  investigated. 

The  Free  Flight  tests,  on  the  other  hand,  lie  on  the 
average,  some  30$  below  the  theory  of  Ref.  25«  (Fig.  4). 

It  may  also  be  noted  that  possibly  10>  of  this  disirepancy 
could  be  due  to  conduction  effects,  as.  discussed  in 
Section  4*3.1.  for  the  test  results  of  Ref. 6.  Bearing 
this  in  mind,  together  with  the  difficulty  of  making  Free 
Flight  heat  transfer  measurements,  the  Free  Flight  data 
agree  satisfactorily  with  the  Wind  Tunnel  results. 

The  #orresponding  results  for  turbulent  flow  over 
hemispheres  and  nose  caps  are  shown  in  Figs  5  and  6,  for 
both  Wind  Tunnel  and  Free  Flight  tests.  In  Fig.  5  both 
the  Nusselt  and  Reynolds  numbers  have  been  based  on  the 
distance  measured  from  the  stagnation  point,  s,  (c.f.Fig  1) 
and  the  test  data  are  compared  with  turbulent  flat-plate 
theory.  Although  the  simple  4/5ths  powor  law  represents 
the  trend  of  the  data  well,  tho  bulk  of  the  results  are 
perhaps  better  correlated  by  a  line  about  30$  below  the 
flat  plate  theory. 

In  Ref.  26  a  simple  formula  for  turbulent  flow  over 
hemispheres  is  presented  and  the  results  are  compared  with 
this  theory  in  Fig.  6.  Here,  in  order  to  account  for  the 
presence  of  the  initial  run  of  laminar  flow,  both  the 
Nusselt  and  Reynolds  numbers  have  been  based  on  the 
distance  measured  from  the  transion  point,  s^(c«f’.  Fig  1) 

It  will  be  observed  immediately  that  the  equation  suggested 

in  Ref.  26  is  not  supported  by  the  experimental  evidence 
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cited  here.  In  fact,  the  line  30%  below  the  flat  plate 
theory  suggested  above,  in  connection  with  Pig  5,  would 
also  best  correlate  the  results  in  terms  of  the 
parameters  presented  in  Pig.  6.  Since  transition  is 
impossible  to  predict  accurately  when  carrying  out  project 
studies,  we  may  recommend  that,  for  turbulent  flow  over 
hemispheres,  the  best  results  for  either  Wind  Tunnel 
or  Free  Plight  tests  will  be  obtained  by  basing  the 
Nusselt  and  Reynolds  numbers  on  s  (c.f.  Pig*  1)  and  using 
the  simple  flat-plate  equation  with  about  a  30%  reduction 
factor. 

4- 1 • 2.  Results  for  conical  skirts  (Figs.  7-13) 

The  variation  of  Nusselt  number  with  Reynolds  number 
for  laminar  flow  over  conical  skirts  is  shown  in  Pigs.  7 
and  8,  and  9  and  10,  respectively  for  Wind  Tunnel  and 
Free  Plight  tests.  For  the  Wind  Tunnel  results  presented 
in  Pig.  7*  it  will  be  seen  that  the  results  lie  between 
the  flat  plate  and  sharp  cone  theories.  This  result ,aiy 
be  anticipated  intuitively,  for  at  very  large  values  of 
s/R  on  a  blunt  cone,  the  effect  of  tip  blunting  should  be 
insignificant  with  the  flow  approaching  fully  developed 
conical  conditions.  Lees,  in  Ref.  27,  also  shows  this 
effect  as  a  limit  to  his  theory  for  large  s/R  values.  The 
existence  of  this  trend  to  conical  flow  with  increasing 
s/R  values  is  discussed  later  in  Section  4.3.2. 

An  alternative  measure  of  the  Nusselt  and  Reynolds 
numbers  is  obtained  by  basing  them  on  the  distance  measured 
from  the  virtual  apex  of  the  cone,  s„  (c.f.  Pig.  1).  The 
results  based  on  this  distance  are  shown  in  Pig.  8,  and  it 
will  be  observed  that,  for  the  results  of  Ref.  tt,  the 
correlation  is  almost  identical  with  that  of  Pig.  7.  The 
results  of  Ref.  11,  however,  become  more  evenly  distributed 
between  the  flat-plate  and  conical  theories.  However, 
comparison  of  Pigs  7  and  8  would  indicate  that  there  is 
really  nothing  to  choose  between  these  alternative  methods 
of  estimating  the  Nusselt  and  Reynolds  numbers,  at  least 
for  the  Wind  Tunnel  test  results. 

The  corresponding  results  for  Free  Plight  tests  are 
shown  in  Pigs  9  and  10.  Two  complete  sets  of  results  are 
shown,  the  open  symbols  corresponding  to  normal  shock 
wave  conditions  and  the  full  symbols  to  conical  shock 
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wave  conditions  (c.f.  Section  3. 2.1.). 

Examination  of  the  data  presented  in  Pig.  9*  indicates 
that  the  two  methods  suggested  for  determining  the  local 
flow  conditions  correlate  the  results  equally  well  between 
the  flap-plate  and  cone  theories.  Consequently,  providing 
that  the  data  is  consistently  analysed  it  would  appear 
not  to  matter  which  assumption  is  made  for  estimating  tho 
local  flow  conditions.  Examination  of  these  results  in 
conjunction  with  the  nose  to  base  radius  ratios  for  each 
test  (tabulated  in  Fig.  2)  was  also  made. 

This  examination  indicated  (for  the  limited  data  of 
three  tests)  that  the  results  are  not  always  improved  by 
changing  from  conical  shock  wave  conditions  at  small 
R/Rb  values  (e.g.  0.138  for  Ref.  17)  to  normal  shock  wave 
conditions  at  large  values  of  R/Rb  (e.g.  0.333  for  Ref . 1 ) 
So  again,  it  may  be  concluded  that  there  is  little  to 
choose  between  the  two  methods  of  estimating  the  flow 
conditions. 

In  Pig.  lO  the  Free  Plight  results  are  shown  in  terms 
of  the  Nussolt  and  Reynolds  numbers  based  on  the  distance 
from  the  virtual  apex  of  the  cone  skirts.  A  significant 
change  occurs  in  the  correlation  when  this  is  done,  as 
can  be  seen  by  comparing  Pigs.  9  and  10,  The  results, 
for  either  set  of  local  flow  conditions,  are  almost  uni¬ 
formly  distributed  about  the  cone  theory  line.  The  low 
points  of  Ref.  17  in  the  Reynolds  number  range  10^  to 
2  x  106  may  be  due  to  conduction  effects.  With  these  few 
points  disregarded,  one  may  conclude  that,  contrary  to  the 
Wind  Tunnel  results,  tho  Free  Plight  laminar  results  on 
conical  skirts  arc  best  correlated  by  sharp  cone  theory 
with  tho  Nussolt  and  Reynolds  numbers  based  on  sa,  the 
distanco  from  the  virtual  apex  of  the  cone  skirt. 

The  Wind  Tunnol  and  Free  Plight  results  for  turbulent 
flow  over  cone  skirts  are  shown  in  Pigs.  11  to  13.  Both 
normal  and  conical  shock  wave  conditions  are  shown.  As 
for  the  laminar  cases,  there  is  little  to  choose  between 
the  two  sets  of  local  flow  conditions,  except  for  the 
results  of  Ref.  4.  which  aie  better  when  based  on 
conical  shock  wave  conditions.  The  results  for  the  model 
of  Ref.  17,  which  has  a  smaller  R/Rq,  value  than  that  of 
Ref.  4,  are  not  improved  by  selecting  conical  shock  wave 
conditions.  As  with  the  laminar  results,  then,  we  find 
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no  consistent  effect  of  the  nose  to  base  radius  ratio  on 
the  selection  of  the  best  local  flow  conditions  to  use. 

Within  the  scatter  of  the  data  there  is  little  to 
choose  bet ween  the  three  correlations  presented  in  Figs. 

11,  12,  and  13  though  the  latter,  with  the  parameters 

based  on  the  distance  from  the  transition  point 

(s-t  c.f.  Fig.  1 )  is  slightly  better  than  the  other  two. 

It  may  be  concluded,  then,  that  for  turbulent  flow  over 
conical  skirts  the  results  are  most  simply  correlated  in 
terms  of  the  Nusselt  and  Reynolds  number  based  on  either 
of  the  lengths  sa  or  s-^  .  Furthermore,  with  the  scatter 
in  the  results  no  clear  choice  between  the  cone  and  flat 
plate  theories  can  be  made. 

4.2.  Variation  of  the  laminar  heat  transfer  parameter  NuD/(P  D^/v  )^ 
with  distance  round  a  sphere  surface  (Fig.  14) ♦ 

The  Variation  of  the  laminar  heat  transfer  parameter  with 
s/R  is  shown  in  Fig.  14  for  laminar  flow  Wind  Tunnel  tests  on 
hemispheres  and  cone  nose  caps.  For  comparison,  Korobkin's  theory 
(Ref.  31)  is  also  shown.  It  can  be  seen  that,  as  predicted  by 
this  theory,  the  heat  transfer  parameter  decreased  with  s/R. 
Although  this  parameter  is  independent  of  Reynolds  number  it  is 
not  independent  of  Mach  number  and  wall  temperature  ratio  effects. 
To  compare  the  results  with  theories  such  as  those  presented 
in  Refs.  5,  27  and  28,  for  example,  would  involve  evaluating 
these  theories  for  each  set  of  conditions  pertaining  to  the 
experiments  under  consideration.  Therefore,  although  this 
parameter  is  frequently  used  to  present  heat  transfer  results 
obtained  from  laminar  flow  tests ,  it  was  abandoned  in  favour  of 
another  parameter,  discussed  below. 

4.3.  Correlations  of  the  ratio  of  Compressible  to  Incompressible 

Nusselt  numbers  with  s/R.  (Figs  15  -  22). 

4.3.I.  Results  for  hemispheres  and  cone  c?.ps  (Figs.  15  -  18) 

The  ratios  of  compressible  (i.e.  experimental)  to 
incompressible  Nusselt  number  plotted  against  s/R  for 
hemispheres  and  nose  caps  are  shown  in  Figs.  1 5  —  18.  As 
shown  in  Appendix  IV^the  parameter  plotted  throughout 
these  Figures  (excluding  Fig.  17)  should  be  independent 
of  all  the  parameters  M-j ,  Tw/T-j,  s/R  and  Res  •  In  Fig.  15 
the  results  for  laminar  flow  wind  tunnel  tests  are  seen 
to  be  evenly  distributed  about  the  theoretical  value  of 
unity  for  this  parameter.  Inspection  of  the  values  of 
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Tw/T-j ,  Res  and  M-j  for  all  the  points  also  showed  tho  results 
to  he  independent  of  these  parameters  although  these  are 
not  noted  in  Fig.  15* 

The  corresponding  Free  Flight  results,  shown  in  Fig.  16, 
as  might  be  anticipated  from  the  results  in  Fig.  4 
(discussed  in  Section  4.2.)  lie  below  the  theoretical  value 
of  unity.  However,  the  data  presented  in  terms  of  this 
parameter  is  again  independent  of  s/R.  The  points  obtained 
from  Ref.  6,  close  to  the  stagnation  point  are  seen 
to  be  rather  low.  As  shown  in  Appendix  II,  if  longitudinal 
conduction  effects  are  ignored  in  estimating  the  Nusselt 
number  this  quantity  may  be  under  estimated.  The 
correction  for  heat  conduction  effects  was  applied  to  the 
results  of  Ref.  6  (at  a  flight  time  of  5  secs)  the  result 
of  which  can  be  seen  in  Fig.  17.  It  can  be  seen  that  the 
correction  to  the  Nusselt  number  can  reach  about  lO^fe 
which  could  improve  the  low  data  points  of  Fig.  16. 

(it  should  be  noted  that  the  data  for  Ref.  6  shown  on 
Fig.  16  corresponds  to  various  flight  times  so  that  tho 
above  statement  is  only  qualitative)* 

The  results  for  turbulent  flow  over  hemispheres  and 
nose  caps  for  both  Wind  Tunnel  and  Free  Flight  tests  are 
shown  in  Fig.  18.  As  with  the  laminar  flow  results  these 
are  again  seen  to  be  essentially  independent  of  s/R  and 
examination  of  the  conditions  for  each  individual  point 
on  the  Figure  showed  the  results  to  be  independent  of 
Mach  number,  Reynolds  number  and  wall  temperature  ratio. 

4.3.2.  Results  for  conical  skirts 

The  Nusselt  number  ratios  obtained  from  Wind  Tunnel 
tests  on  conical  skirts  with  laminar  flow  are  shown  in 
Fig.  19.  For  the  results  of  Ref.  10  (oircle  points)  there 
•  is  a  distinct  trend  for  the  results  to  approach  the 
theoretical  conical  value  as  s/R  becomes  large.  The 
results  for  Ref.  11  (triangle  points),  however,  do  not 
exhibit  this  trend.  If  the  values  of  R/r^  for  these  two 
tests  are  compared  (see  Fig.  2)  it  will  be  observed  that 
the  model  of  Ref.  11  is  effectively  twice  as  blunt  as 
that  tested  in  Ref.  10.  Hence,  the  trend  towards  conical 
flow  conditions  would  appear  to  be  determined  primarily 
by  the  nose  to  base  radius  ratio  in  addition  to  the  value 
of  s/R. 
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The  corresponding  Free  Flight  results  are  shown  for 

normal  and  conical  shock  wave  conditions  respectively  in 

Figs.  20  and  21.  Comparison  of  the  two  figures  supports  the 

earlier  conclusion  (made  in  Section  4* 1.2.)  that  basing  the 

local  flow  conditions  on  either  conical  or  normal  shock  wave 

conditions  makes  little  difference  to  the  correlations.  As 

found  in  the  Nusselt  number  -  Reynolds  number  plot  (Fig.  9)> 

the  results  fall  between  the  cone  and  flat  plate  theories 

except  at  .the  lower  values  of  s/R  (  <11).  Here  the  effects 

of  nose  blunting  may  bo  thought  to  be  responsible  for  the 

spread  in  the  data.  However,  the  results  of  Ref.  1  (which 

correspond  to  r/R'  *=  0.333)  lie  consistently  lower  than 
Id 

those  of  Ref.  4>  which  correspond  to  R/Rt  “  0.235*  Also, 
the  points  obtained  from  Ref.  17  (R/Rb  *=  0.138)  lie  slightly 
lower  than  those  of  Ref.  1  so  that  there  is  no  consistent 
trend  of  the  data  in  this  region  of  s/R  with  nose  blunting. 
The  lower  points  (Refs.  17  end  1)  may  be  low  due  to 
conduction  effects  previously  described,  or  they  may  be  due 
to  errors  in  the  pressure  distribution  (c.f.  Section  3*1*) 

It  can  also  be  observed  that,  unlike  the  Wind  Tunnel 
results  (c.f.  Fig.  19),  no  trend  of  the  Nusselt  number  ratio 
with  s/R  is  present.  Examination  of  the  Mach  number, 
Reynolds  number  and  wall  temperature  ratios  for  the  points 
presented  on  Figs.  20  and  21  again  show  the  Nusselt  number 
ratio  parameter  to  be  independent  of  these  variables. 

The  corresponding  Wind  Tunnel  and  Free  Flight  tests 
for  turbulent  flow  are  shown  in  Fig.  22.  Here  again  little 
difference  is  seen  between  the  two  methods  of  estimating 
the  local  flow  conditions.  Also  as  found  for  the  Free 
Flight  laminar  flow  results  there  is  a  large  scatter  in 
the  data  for  values  of  s/R  <12.  Because  of  the  small 
difference  between  cone  and  flat  plate  theories  (about  15$) 
for  turbulent  flow,  the  scatter  in  the  data  does  not  allow 
the  detection  of  a  trend  towards  the  conical  value  of  the 
Nusselt  number  ratio  for  large  s/R  values,  as  was  possible 
for  the  laminar  flow  Wind  Tunnel  results.  However, 
examination  of  the  Mach  number,  Reynolds  number  and  wall 
temperature  ratios  for  these  points  show  the  parameter 
to  be  essentially  independent  of  these. 
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4. 4*  Reynolds  Analogy  Factor  Results  (Figs.  23  -  33) 

4.4.I.  Introductory  Discussion 

The  details  of  the  methods  used  for  estimating 
Reynolds  Analogy  Factor  are  described  in  Appendix  V.  For 
hemispheres  the  first  method  used  was  to  base  the  local 
skin  friction  coefficients  on  simple  formulae  of  the  type 
used  in  the  Nusselt  numb er-Reyno Ids  number  correlations  of 
Figs. 3  &5,In  the  case  of  the  hemispheres,  only  the  laminar 
flow  Wind  Tunnel  tests  were  examined  for  evidence 
supporting  Reynolds  Analogy. 

The  second  method  was  to  derive  an  expression  relating 
the  local  skin  friction  coefficient  cn  an  axially  symmetric 
body  with  either  laminar  or  turbulent  flow  to  that  on  a 
flat  plate  in  terms  of  a  pressure  gradient  parameter,  \)r  . 
The  details  of  this  analysis  are  outlined  in  Appendix  VI. 
This  method  of  evaluating  the  skin  friction  coefficients 
was  then  used  as  a  basis  for  estimating  Reynolds  Analogy 
Factor  as  described  in  Appendix  V.  This  second  method, 
owing  to  the  length  of  time  involved  in  its  calculation, 
was  only  applied  to  the  laminar  flow  results  as  a  sample 
test  of  the  method. 

Many  authors  (e.g.  Lees  (Ref.  27)  ),  Vaglio-Laurin 
(Refs.  28,  29) ,  and  others  reviewed  by  Phillips  (  Ref. 3*) 
have  derived  similar  pressure  gradient  parameters  to  the 
one  derived  in  the  present  report,  but  their  analyses  are 
restricted  to  the  hypersonic  'cool  wall'  approximation. 
However,  the  results  presented  here  are  not  represented  by 
this  approximation,  and  thus  a  special  analysis  was 
carried  out  to  determine  i|r  .  It  should  be  noted  that  the 
method  of  Steine  and  Wanlass  (Ref.  5)>  though  not 
restricted  to  'cool  walls'  is  laborious  to  evaluate,  and 
was  -ihorefore  unsuitable  for  use  in  the  present 
investigation  . 

For  estimating  Reynolds  Analogy  Factor  on  cone  skirts, 
the  local  skin  friction  coefficients  were  based  throughout 
on  the  sharp  cone  equations. 
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4*4*2.  Reynolds  Analogy  Factor  for  hemispheres  and  nose  caps 

(Figs.  24-27) 

Reynolds  Analogy  Factors  for  laminar  flow  over 
hemispheres  as  found  using  the  two  methods  discussed  in  the 
previous  section  are  shown  in  Figures  24  to  27  for  the  Wind 
Tunnel  tests.  Figs  24  and  25  show  the  results  plotted 
against  local  Mach  number,  whereas,  in  Figs.  26  and  27  they 
are  plotted  against  local  Reynolds  number.  The  results  in 
Figs*  24  and  26  show  immediately  the  great  improvement  in 
the  correlation  obtained  by  using  the  pressure  gradient 
parameter  i|r .  However,  if  the  results  based  on  the  use  of 
the  parameter  are  examined  on  a  larger  scale  (Figs.  25 
and  27)  they  are  seen  to  lie  almost  wholly  below  the  value 
of  Pr"^/suggested  by  Colburn  in  Ref.  24.  For  either  method 
used  to  estimate  Reynolds  Analogy  Factors  on  hemispheres 
for  laminar  flow,  the  results  are  seen  to  be  virtually 
independent  of  local  Mach  number  and  Reynolds  number. 
Although  the  use  of  the  pressure  gradient  parameter, +  in 
estimating  Reynolds  Analogy  Factor  brings  the  experimental 
results  more  into  line  with  Colburn's  value  (s  =  1.246)^ 
it  does  increase  the  scatter  among  the  data  points. 

(c.f.  Figs  24  and  25,  for  example). 

4»4»3.  Reynolds  Analogy  Factor  for  conical  skirts  (Figs  28-33) 

Figs  28  and  29  show  the  values  of  Reynolds  Analogy 
Factor  obtained  from  laminar  flow  Wind  Tunnel  tests  on  cone 
skirts.  They  are  plotted,  respectively  against  local  Mach 
number  and  local  Reynolds  number  in  the  two  figures.  The 
data  here  is  so  limited  that  no  trend  of  Reynolds  Analogy 
Factor  with  either  Reynolds  number  or  Mach  number  can  be 
observed.  It  is  worth  noting,  however,  that  the  use  of  the 
pressure  gradient  parameter,^  ,  in  place  of  the  pointed 
cone  factor,  V3>  in  determining  Reynolds  Analogy  Factor 
for  these  tests  might  correlate  the  results  better  with  the 
line  Pr  -3=  1.246.  (This  is  providing  that  the  value  of 
on  cone  skirts  does  not  vary  significantly  from  that  shown 
in  Fig.  23). 

The  corresponding  Free  Flight  results  for  laminar  flow 
over  cone  skirts  are  shown  in  Figs.  30  and  31.  The  results 
based  on  both  methods  of  estimating  the  local  flow 
conditions  are  shown.  Although  in  both  Figures  the  data  is 
moved  to  the  right  in  going  from  normal  to  conical  shock 
wave  conditions,  the  value  of  Reynolds  Analogy  Factor  for 
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each  point  remains  virtually  unaltered.  The  data  is  seen 
to  fall  below  the  theoretical  value.  This  is  partly  due 
to  the  low  values  of  the  heat  transfer  data  observed  in 
Pig.  9  and  also,  as  suggested  for  the  Wind  Tunnel  tests, 
the  results  might  be  improved  by  the  introduction  of  the 
pressure  gradient  parameter  in  the  place  of  the  V" 3  factor. 
Again  the  data  shows  no  definite  trend  with  either  local 
Mach  number  or  local  Reynolds  number. 

To  conclude  the  Reynolds  Analogy  Factor  investigations, 
the  results  for  turbulent  flow  over  cone  skirts  are  shown 
in  Figs.  32  and  33,  plotted  against  local  Mach  number 
and  Reynolds  number  respectively.  Again  the  results 
obtained  are  independent  of  the  method  of  estimating  the 
local  flow  conditions  and  also  of  the  Reynolds  and  Mach 
numbers.  This  latter  result  was  observed  for  the  Free 

Flight  tests  on  pointed  cones  investigated  in  Ref. 21. 
Estimation  of  Wall  Temperatures  in  Free  Flight  for  Thin  Skin  Models 

Theoretical  estimates  were  made  of  Hie  wall  temperatures 
attained  during  the  flight  of  four  rocket  propelled  blunt  cone 
models  and  one  hemisphere  model.  These  were  based  on  the  method 
of  Hill  (Ref.  32)  with  the  heat-transfer  coefficients  determined 
from  the  semi-empirical  equations  duo  to  Monaghan  (Ref.  33) 

(Details  of  the  calculations  are  presented  in.  Appendix  VII. 

The  ratio  of  the  measured  to  theoretical  wall  temperatures 
for  selected  measuring  stations  on  the  models  of  Refs.  1,2,  4 
and  6  are  plotted  in  Fig.  34  against  the  reduced  flight  time  t 
(defined  on  the  Figure). 

Considering  the  simplicity  of  Monaghan's  equations 
(Ref.  33)  the  results  are  very  good,  almost  all  of  them  lying 
within  10$  of  the  perfect  correlation;  (i.e.  the  line  1.0  on 
Fig.  33)  The  method  used  here  is  thus  to  be  recommended  for 
estimating  wall  temperatures  for  thermally  thin  skins. 
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5.  CONCLUSIONS 

5.1 .  Nusselt  number  -  Reynolds  number  Relations 

i)  For  laminar  flow  over  hemispheres  and  nose  caps,  for 

2  z  10^  <  H-*  <  3  x  10^  ,  the  Wind  Tunnel  results  correlate 
s 

well  with  Sibulkin's  formula  expressed  in  the  form 
Nus*  Pr  *  3  „  0.763  Res*  * 

The  corresponding  Free  Flight  results  follow  the  trend 
indicated  by  the  above  relation  but  lie  20  to  25fc  below  the 
values  given  by  this  equation.  The  Reynolds  number  range 
for  these  results  is  1.8  x  104^Res«  6.2  x  10^ 

ii)  For  turbulent  flow  over  hemispheres  and  nose  caps,  both  the 
Wind  Tunnel  and  Free  Flight  data  correlate  best  with  the 
flat  plate  equation  with  the  constant  reduced  by  about 

30^  i.o.  _± 

Nus  Pr  *  (0.7)  x  (0.#296)  Res 

These  results,  when  plotted  in  terms  of  the  parameters  based 
on  the  distance  measured  from  the  transition  point,  are  not 

correlated  by  van  Dries t's  formula,  namely, 

*  *_i  *  4/5 

Nut  Pr  “  0,042  Rot 

iii)  The  results  obtained  from  Wind  Tunnel  tests  on  cone  skirts 
with  laminar  flow  when  correlated  in  terms  of  the  parameter 
Nu*  Pr*  -3  vs  Re  lie  between  the  flat  plate  and  sharp 
cone  theories.  This  result  is  independent  of  the  lengths 
on  which  Nusselt  and  Reynolds  numbers  are  based. 

iv)  The  results  obtained  from  Free  Flight  tests  on  cone  skirts, 

for  laminar  flow,  are  best  correlated  in  terms  of  the  Nusselt 
and  Reynolds  numbers  based  on  the  distance  measured  from 
the  virtual  apex  of  the  cone  skirt  by  the  equation. 

Nu*a  Pr*  «  V3  x  0.332  Rea*^ 
for  106<  Roa  <  2  x  107 

v)  For  both  Wind  Tunnel  and  Free  Flight  tests  on  cone  skirts 

with  turbulent  flow,  tho  scatter  in  the  measurements  is  such 
that  the  data  correlates  equally  well  with  either  the  flat 
plate  or  sharp  cone  theories.  This  is  true  regardless  of  the 
distance  on  which  the  Nusselt  and  Reynolds  numbers  are 
based. 
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vi)  For  the  Free  Flight  tests  on  conical  skirts,  the  correlations 
of  Nusselt  number  with  Reynolds  number  are  independent  of 
whether  the  local  flow  parameters  are  based  on  normal  or 
conical  shock  wave  stagnation  conditions.  This  is, true  for 
both  the  laminar  and  turbulent  flow  results. 

Results  in  terms  of  the  Compressible  to  Incompressible  Nusselt 

number  Ratios 

i)  For  both  laminar  and  turbulent  flow  over  hemispheres  and 
nose  caps,  the  experimental  results,  when  reduced  to  the 
form  (Nug/Nui)  (T*/T-|)a  aro  independent  of  s/R.  This  is 
so  when  Nuj_  is  based  on  Sibulkin’s  equations  for  laminar 
flow,  or  with  Nu.  based  on  the  flat  plate  relation  for 
turbulent  flow.  Furthermore,  with  a  =  0.13  and  0.64  for 
laminar  and  turbulent  flow  respectively,  the  parameter 
(Nus/Nuj_)  (T*/Ti)a  is  essentially  independent  of  the  local 
values  of  wall  temperature  ratio,  Mach  number  and  Reynolds 
number. 

ii)  For  cone  skirts  the  results  obtained  from  laminar  flow 
Wind  Tunnel  tests,  expressed  in  terms  of  the  parameter 
(Nug/Nui)  (TVT^O-13  aro  independent  of  the  local  values  of 
wall  temperature  ratio,  Mach  number  and  Reynolds  number. 
Furthermore  the  trend  of  this  parameter,  when  based  on 
flat  plate  theory,  towards  tho  sharp  cone  value  of 
for  large  values  of  s/R  appears  to  depend  primarily  on  the 
nose  to  base  radius  ratio.  Conclusive  evidence  of  this 
effect  could  not  be  found  because  of  The  limited  number  of 
tosts  investigated. 

iii)  For  conical  skirts  the  results  obtained  from  laminar  flow 

Free  Flight  tests,  when  expressed  in  terms  of  the  parameter 
(Nus/Nuj_)  (T*/T-|)0#64  show  no  distinct  trend  towards  the 
sharp  cone  limit  at  large  values  of  s/R.  Again  this  parameter 
is  essentially  independent  of  the  local  values  of  wall 
temperature  ratio,  Mach  number  and  Reynolds  number. 

iv)  For  the  turbulent  Free  Flight  and  Wind  Tunnel  test  results 
on  cone  skirts,  expressed  in  the  form  (Nug/Nu^)  (T*/T-|)0.64 
the  scatter  in  the  data  allows  no  trend  between  the  flat 
plate  and  sharp  cone  limits  to  be  observed.  This  parameter 
is  again  found  to  be  independent  of  the  local  values  of  wall 
temperature  ratio,  Mach  number  and  Reynolds  number. 
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v)  For  both  laminar  and  turbulent  flow  Free  Flight  tests 
conclusions  5  (ii)  to  (iv)  are  independent  of  whether 
the  data  is  based  on  normal  shock  wave  or  conical  shock 
wave  stagnation  conditions. 

vi)  For  Free  Flight  tests  on  cone  skirts  the  parameter 

(Nus/Nuj_)  (T*/T-|)a  has  a  large  scatter  for  values  of  s/R 
less  than  twelve.  This  scatter  could  not  be  correlated 
with  either  the  nose  to  base  radius  ratios  nor  was  it 
completely  associated  with  conduction  effects  along 
the  model  surfaces*  Also  some  of  this  scatter  may  be  due 
to  errors  in  estimating  the  pressure  distribution,  as 
this  is  not  always  quoted  precisely  in  Free  Flight  test 
reports. 

5.3.  Reynolds  Analogy  Factors 

i)  For  the  data  investigated  Reynolds  Analogy  Factor  is 

virtually  independent  of  both  local  Mach  number  and  local 
Reynolds  number. 

il)  For  hemispheres  with  laminar  boundary  layers,  Reynolds 

Analogy  Factor  is  best  determined  from  experimental  heat 
transfer  measurements  when  the  local  skin  friction 
coefficients  are  corrected  for  pressure  gradient  and  body 
shape  effects. 

5.4.  Wall  temperatures  on  Thin  Skin  Free  Flight  Models 

For  the  cases  investigated,  9 0$  of  the  measured  wall 
temperatures  are  within  10$  of  the  theoretical  estimates  for 
hemispherical  bodies  and  cone  skirts. 
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APPENDIX  I 

Estimation  of  probable  »^rors  in  n  iasured  Stanton 
numbers  in  Free  Flight  tests  and  the  corresponding  errors 

in  Nusselt  numbers 


In  Appendix  3  of  Ref.  21,  it  is  shown  that  in  the  absence 
of  errors  in  local  pressure  measurements,  the  probable  error  in 
the  measured  local  Stanton  number  is  given  by, 


t  otal 


P 401  +  0.16  ftfco  +(T4(l-r)  0.01)3  , 

to  Y  +  0  00132  1 

L  (Tr  -  ^ 

v1TJ  3  J 

++ 


I  (1) 


(The  errors  assumed  in  individual  quantities  are  tabulated 
at  the  end  of  this  Appendix) • 


Since,  in  the  present  report,  we  deal  exclusively  7/ith  Nusselt 
numbers  we  now  derive  an  expression  relating  errors  in  Nusselt  number 
to  those  in  Stanton  number. 


Noting  that, 

Nu  Pr.  Re. 

it  is  easily  shown,  using  the  relations  in  Appendix  3  of  Ref«  21, 
that  the  probable  error  in  the  Nusselt  number  is  given  by, 


I  (2) 


thus,  ignoring  errors  in  the  local  flow  conditions  implied  by 
pressure  distribution  errors,  we  find, 


/ crNu 
^  Nu 


total 


fa ST 
V  ST  A  otal 


•  I  (3) 


so  that  equation  I  (l)  also  gives  the  probable  error  in  the 
Nusselt  number. 

It  is  also  necessary  to  show  what  the  probable  errors  in  Nu  are, 
since  this  parameter  is  also  used  extensively  in  the  correlations. 


We  have. 


Nu  =  Nu 


(ki  /k i‘ 


.1  (4) 


and,  using  the  relations 

k  ~  T'1'  j  p  ~  l/T  (Since  p  =  const, .across  the  boundary  layer) 


++  Temperatures  in  eqn.  I  (l)  are  in  °R, 
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equation  I  (4)becomes : 

Nu*  =  Nu  (p4/p*  J1 


I  (5) 


Consequently,  using  the  previously  quoted  method,  we  find 


I  (6) 


Now,  inspection  of  the  flow  data  for  the  reports  investigated 
here,  shows  that. 


i 2  r<r  pi/  p  *  "i  roifi 

J <<:  LNuJ 


I  (7) 


Consequently, 

* 

0£ 


(<f  Nu  *  \  M  fa-  Nu  \ 

\  Nu  *  /  V  Nu/ 


I  (8) 


•  3E 

Finally,  the  errors  in  Re  and  Pr  due  to  errors  in  measured 
wall  temperature  can  he  shown  to  have  no  significant  effect  on  the 
correlations. 


Table  of  probable  errors  (After  Ref.  22) , 


Quantity 

Error 

T 

w 

20°  R 

T  shield 

20°  R 

Poo 

0.013  p 

1  OO 

T 

1°R 

00 

u 

oo 

4  ft/ sec 

dT  /&T 
w' 

1°  R/sec 

X 

III 

0.03  x 

P1/  poo,  Ti / T  eo,  Ui j Uoo 

0 

€s 

Hrm — ~  /a  ruv - 

0.02 
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Appendix  II 

Estimation  of  bho  effect  of  longitudinal  conduction 

on  Nusselt  number. 


In  the  absence  of  circumferential  heat  flow  and  radiation 
the  heat  balance  equation  for  a  typical  surface  element  may  be 
written 


\  -  (p  °  iT„  -  ks  -1  fi-  p  »Tw  .  a2Tw  ] 

3T  1/  38  —  +  37“  J 


II  (1) 


where  q^  =  convective  heat  transfer  rate  per  unit  area 
per  unit  time# 

y  =  body  perimeter  in  planes  normal  to  the  axis. 

Defining  a  local  Nusselt  number  by, 

Nu±=  q  s/kt  (T  -  T  )  . II  (2) 

^a  '  v  r  w  7 

equation  II  (l)  may  be  written  as 

Nu±=  (Nu)  -  k  'V  s  ~  tt 

7  exp  s  Q  . II  (JJ 

kTJr  -  t  ) 

where  Q  is  the  term  in  curly  brackets  in  equation  II  (l)  and  (Nu)  ^ 
is  the  experimentally  determined  Nusselt  number  neglecting  conduction 
and  radiation  effects,  i.e., 

(Nu)exp  =  (Pct)s*  z  fT  T'  ~y  . II  (4) 

If  the  Nusselt  number  is  based  on  enthalpy  in  place  of 
temperature,  equation  II  (3)  becomes 

Nu±  =  (Nu)  -kvsCQ  \ 

exp  s  p  ,  . II  (5; 

ki (i  -  i  ) 

x  r  w' 


and  C  is  such  that 
P 

C  (T  -  T  )  =  (i  -  i  )  . II  (6) 

p  N  r  w7  r  w7  •  \  / 

Consequently,  depending  upon  the  sign  of  Q  in  equation  II  (5) 
the  experimentally  determined  Nusselt  number  can  either  over  - 

V 

or  underestimate  the  true  Nusselt  number,  if  conduction  is  ignored. 
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Appendix  III 

Summary  of  the  equations  used  in  the  Nusselt  number  - 

Reynolds  number  correlations. 


1 .  Hemispheres  and  nose  caps. 

For  laminar  flow  over  hemispheres  an  equation  relating  the  local 
Stanton  and  Reynolds  number  is  given  in  Ref.  25  in  the  form, 


ST  =  0.763  Pr  "°*6  Ros 


III  (1) 


Replacing  ?r  0,6  by  Pr“  2^3 


and  introducing  intermediate 
enthalpy  conditions  to  account  for  compressibility  effects,  this  equation 
becomes 


* 

S_  P 
T  r 


*  2/3 


0.763  Re 


x 

*  -  2 


Noting  that  the  Stanton  and  Nusselt  numbers  arc  related  "by 


*  *  *  * 

N  =  Sm  P  Re 
u  T  r  s 

s 


m  (2) 


the  Nusselt  number  for  hemispheres  in  compressible  flow  is  found  to  be 

M/3 


* 

N  P 
u  r 
s 


0.763  Ro  *  ^ 

S 


in  (3) 


For  turbulent  flow  over  hemispheres,  the  equation  given  in  Ref, 26, 
in  terms  of  local  flow  conditions  is 


Sm  »  0.042  P  2/3  Re  “  1//5 

T  r  s 


hi  (4) 


Following  the  same  procedure  as  for  laminar  flow,  the  Nusselt 
number  based  on  intermediate  enthalpy  conditions  for  turbulent  flow  over 
hemispheres  is  found  from  III  (4)  to  be 


*  *  -  1/3 


N  Pr 
u 

s 


0.042.  Re 


„  *4/5 


in  (5) 


To  correct  for  the  initial  run  of  laminar  flow  the  Nusselt  and 
Reynolds  numbers  arc  based  on  the  distance  measured  from  the  transition 
point  (Fig.  1 )  so  that  equation  III  (5)  becomes 


# 

N  P 
ut  r 


*  -  1/3 


=  0.042  R 


*  4/5 


III  (6) 


*  _  1/3 

++  This  has  been  done  merely  to  preserve  the  parameter  P^  '  J  throughout 
the  correlations. 
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2.  Flat  Plate 

For  laminar  flow  the  Nusselt  and  Reynolds  numbers,  based  on 
intermediate  enthalpy  conditions  are  related  by 


N 


u 


P 


r 


III  (7) 


S 


and  for  turbulent  flow  they  arc  related  by 


'u 


* 


IH  (8) 


with  Nu  and  Re  based  on  either  s  or  s^  respectively. 

3.  Cone  Skirts 

For  cone  skirts  the  flat  plate  equations  are  used,  replacing 
the  constants  by 

•  <  /  < 

V3  x  0*332  for  laminar  flow 
and  1,15  x  0.0296  for  turbulent  flow. 


In  addition  to  basing  the  Nusselt  and  Reynolds  numbers  on  the 
distances  sand  (for  turbulent  flow),  they  may  also  be  based  on  the 
distance  measured  from  the  virtual  apex  of  the  cone  skirts,  sa.  (See 
Fig.  1.) 
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Appendix  IV 

Relations  “between  the  compressible  and  incompressible 

Nusselt  numbers  for  laminar  and  turbulent  flew. 


For  a  flat-plate,  in  either  laminar  or  turbulent  flow  the  local 
Nussolt  number  and  local  Reynolds  number  are  related  in  incompressible  flow 
by  an  equation  of  the  form, 

Nu.  ~  Pr  ^3  (Re  )n  IV  (l) 

1  1  1 

with  n  =  4/5  for  turbulent  flow  and  g-  for  laminar  flow. 


Following  tho  reference  temperature  correction  technique  due  to 
Eckert  (Ref.  20),  we  first  find  an  'intermediate*  Nusselt  number,  given 
by  a  relation  of  tho  form, 

Nu  *  ~  Pr  *  1//3  (Re  *)n  _ IV  (2) 

* 

Where  the  flow  properties  are  evaluated  at  a  temperature  T 

* 

corresponding  to  the  enthalpy  i  given  hy 

i*  =  i  +  0.5  (i  -  i  )  +  0,22  (i  -  i  )  _  IV  (3) 

1  w  1  r  1  — 

* 

The  Nusselt  number  in  compressible  flow  is  related  to  Nu  by, 

(Nu)  c  -  Nu  (k*/k4  )  _  IV  (4) 

Hence,  combining  equations  IV  (l),  (2)  and  (4),  the  ratio  of 
compressible  to  incompressible  Nusselt  numbers  is  found  qs, 

Nuc/Nui  =  (Pr*/Prt  )^3  (Re*/^.  )n  ^  \  ^  - IV  ^ 

finally,  making  use  of  the  following  simplifications, 


Pr 


Pr. 


T 


0.7 


nO.723 


then  equation  IV  (5)  becomes, 


(t7t4  ) a 


IV  (6) 


where 


a  =  0.13  for  laminar  flow 

a  =  O.64  for  turbulent  flow. 
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Now,  if  for  hemispheres  in  laminar  flow,  the  experimental  values  of 
* 

Nug  correlated  perfectly  with  the  equations- 

Nu  *  =  0.763Pr  *  ^3  Re  *  *  _ IV  (7) 

S  o 

then,  providing  Ihn  is  determined  from  the  incompressible  form  of 
the  above  equation,  the  experimental  results  plotted  in  the  form 
(ifu  /Nu.)  (T*/t.,)  0,13  vs.  s/R  should  lie  about  the  line  1.0. 

Similarly,  if  for  hemispheres  in  turbulent  flow,  the  flat'  plate 

*  *  -  l/3 

equation  best  correlates  the  results  in  terms  of  Nug  Pr  '  vs. 

Rog  5  then  the  parameter  (Nug/Nu^)  (T  /Tt)  ^  should  collapse  the 
experimental  data  about  the  line  unity,  when  plotted  against  s/R. 

This  is  providing,  of  course,  that  Nu^  is  based  on  the  flat  plate 
equation 

Nu.  =  0.0296  P  *  1//3  Roq4/5  _ IV  (8) 

i  r  •  s 

For  conical  skirts  one  may  expect,  for  large  values  of  s/R  and 
slight  degrees  of  nose-blunting  (i.o.  r/Rj.<<1),  that  the  heat  transfer 
conditions  will  approach  those  for  a  pointed  cone.  In  the  correlations 
of  the  parameter  (Nu  /Nu. )  (T*/T-j)a  vs.  s/R  for  cone  skirts,  the  values 

of  Nu^  wore  determined  from  flat  plate  equations.  Consequently,  in 
regions  where  the  f 1 ow  approaches  pointed  cone  conditions  the  parameter 
(Nu  /Nu.)  (T*/Ti )a  would  take  a  value  different  from  unity.  Intuitively 
this  parameter,  for,  cone  skirts,  should  be  in  the  ranges 

1.0  <  (Nu  /Nu.)  (T*/Ti)  0*13  *  *0  for  laminar  flow 

S  X 

and  1,0  <  (Nu  /W )  (T*/ti)  <  1 .15  for  turbulent  flow. 

Finally  we  note  that,  by  using  the 
technique,  the  function  (T  /t^ ) a  should 
and  wall  temperature  effects. 


intermediate  enthalpy  correction 
account  for  both  Mach  number 
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APPENDIX  V 

Derivation  Qf  Reynolds  /analogy  Factor  from 
pleasured  values  of  gtanton  number 


1 )  Laminar  Flow 

Reynolds  Analogy  Factor,  S,  is  related  to  the  Stanton  number 
and  skin  friction  coefficient  by  the  equation 

S  =  ST/(Cf/2)  _ V  (1) 

Noting  that  Cf,  for  compressible  flow,  can  be  expressed  in  the 

forms- 

Cf/2  "  (P*  /P,)  0  (Re*)"*  _ V  (2) 

Reynolds  Analogy  Factor  becomes 
s  -  ST  /  [(P*  /P^  )  c  (Rqs*  ]  _ V  (3) 

where  ST  is  the  experimentally  determined  Stanton  numbe?:  and  c  takes 
the  values  tabulated  at  the  end  of  this  Appendix. 


2)  Turbulent  Flow 


Assuming  the  Prandtl-Schlichting  formula  for  the  compressible 
skin  friction  coefficient  in  turbulent  flow,  (Ref  37),  equation  V  (1) 
gives  Reynolds  Analogy ’Fan tor  in  the  form 


S  -  Sn 


/  £(p  /  Pi  )  o  (l°gl0 


-2.45  i! 

J 


v  (4) 


whero  c  takes  the  values  tabulated  below. 


3)  Values  of  the  constant,  c 


FLAT  PLiTE 

SHARP  CONE 

HEMISPHERE  j 

Laminar  Flow 

0.332 

/  3  x  0.332 

(0.332)  /if  \ 

Turbulent  Flow 

0.144 

1.15  x  0.144 

(0.144)  /* 

where  is  the  pressure  gradient  parameter  derived  in 

Appendix  VI 
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For  hemispheres  with  laminar  flow  both  the  'flat-plate'  and 
'hemisphere'  values  of  c  were  used  in  the  Reynolds  Analogy  Factor  calculations. 

For  hemispheres  with  turbulent  flow  the  1  flat-plate '  values  of 
c  were  used  and  for  the  cone  skirts,  in  all  cases  the  'sharp  cone'  values 
of  c  were  used  for  determining  Reynolds  Analogy  Factor. 
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Appendix  VI 


The  Effect  of  Pressure  Gradients  on  Skin  Friction  Coefficients 

for  Laminar  and  Turbulent  Flow  on  Blunted  Cones. 


The  skin  friction  coefficient  for  turbulent  flow. 

It  is  shown  in  Ref.  34,  that  by  using  a  modified  form  of 
the  Stewartson  -  Illingworth  transformation  introduced  by  Spence 
in  Ref.  35  »  the  momentum  equation  for  axially  symmetric  turbulent 
flow  may  bo  expressed  in  the  form 8- 


d9  + 
ds 


dU  +  1  dr  ] 
ds  r  ds  j 


where 


c  (1  +  m)  (£)a.  (g)  P 
_ VI  (1) 


0  -  0  <u1i7r0)m 

e  -  G  (s)k 

u1  =  u1  («r 1/2 

Q  o  momentum  thickness 

B  «  (t/V  Hi  +  "1 

EL  n  shape  parameter  for  incompressible  flow 


s  ■  VTo 

5  "  TX 

c,m,a,  (3  =  constants 

(The  other  symbols  are  as  defined  for  the  main  report). 


For  the  isothermal  wall  case,  in  which  we  arc  interested, 
Spence  showed  that  k  and  1  take  the  values 


k  =  1  -  r/2  +  l/(y-  1) 

1  -  1  _  VI  (2) 


Equation  VI  ( 1 )  may  be  integrated  to  give 
0 

where 
and 

Sow,  it  is  also  shown  in  Ref.  35  that  the  local  skin  friction 
coefficient  may  be  found  from  the  relation 


c(  1  +  m)  [  ?(3)°  _  VI  (3) 

f  U)  1c 

f(s)  -  (m/.  fc)6  (z)~P  (r)m  +  1 

6  =  a+  (l  —  l)  B/2 
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cf  =  c.e 

2~  ©  _  VI  (4) 

where  @  »  (£)a(g)  9  _ VI  (5) 

Finally,  combining  equations  VI  (3),  (4)  2nd.  (5),  the  local 
skin  friction  coefficient  is  found  to  be 


where 

Z  (s) 


,\c 


(1  + 

(l+m) 


m 


V< 


.  1  /ra+1  VI 


2«  The  skin  friction  coefficient  for  laminar  flow. 


(6) 

(7) 


A  completely  analogous  method  to  that  outlined  above  for  finding 
the  turbulent  skin  friction  coefficient  may  bo  applied  to  the  laminar 
flow  case.  In  order  to  do  this  it  is  only  necessary  to  adjust  the 
values  of  the  constants  a  ,  p  ,  B  ,  m  and  c  and  to  make  a  simple 
postulation  discussed  below.  It  should  be  not  2d  that  B  ,  a  end  p 
depend  upon  the  value  of  m  (see  Ref. 35)  and  that  c  and  m  are  the 
constants  in  the  equation  relating  the  friction  coefficient  on  a 
flat-plate  to  the  momentum  thickness  Reynolds  number,  i.e. 


ff  =  c  (Re0  )"m 
2 


VI  (8) 


For  the  turbulent  flow  analysis,  Spence  showed,  by  assuming 
a  quadratic  temperature  -  velocity  relationship,  which  is  not 
restricted  tc  unit  Prandtl  number,  that  the  compressible  and 
incompressible  shape  parameters  arc  related  by 


H  =  (Tw/T1 )  H±  +  (Tr/T1 )  -  1  _  VI  (9) 

Now,  for  laminar  flow  with  unit  Prandtl  number,  it  can  bo 
shown  (e.g.  Ref.  36)  that  the  temperature  -  velocity  relationship 
across  the  boundary  layer  is  a  quadratic  of  the  fom 


T  -  TW  +  (’o -*.)(§>  <*0-h)(2 


a, 

) 


VI  (10) 
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This  becomes  identical  with  the  relation  assumed  for  turbulent 
boundary  layers  which  leads  to  equation  VT  (9)  if  Tq  is  replaced 
by  Tr. 


Thus,  if  wo  make  this  modification  to  equation  VI  (10) 
without  attempting  a  rigorous  justification,  the  analysis  for 
the  laminar  boundary  layer  is  formally  identical  with  that 
outlined  for  the  turbulent  boundary  layer.  Hence,  equation  VI 
(6)  gives  the  skin  friction  coefficient  for  either  laminar  or 
turbulent  flow  depending  upon  the  values  sclectod  for  the  constants. 

I 

Derivation  of  the  pressure  gradient  parameter, ,  relatin':  skin 

friction  coefficients  on  an  axially  symmetric  body  to  those  on  a  flat 

plate. 


The  function  z  (s)  defined  in  equation  VII  (7)  may  be 
written  in  the  form 


■  to  -  a  (p'/p , )  (pv//rnA.+i  (s)"/^1 


where  a  _ 


f.  C  (1+m)  1//m+1  j  (m+1) 


VI  (11) 
VI  (12) 


We  note  that  with  m  *=  -J-  and  c  =  0.0128  (turbulent  flow) 

A  takes  the  value  0*0296  and  that  with  m  =  -g-  and  c  =  0.2205 
(laminar  flow)  A  takes  the  value  0.332.  Thus z  (s)  may  be 
written  in  terms  of  the  compressible  skin  friction  coefficient 
fer  a  flat  plate  as  follows, 


(s) 


■(.*) 


A-m/m+1 


VI  (13) 
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If  the  value  of  z  (s)  found  in  equation  VI  (13)  is  substituted  into 
equation  VI  (6),  the  relationship  between  the  skin  friction  coefficients 
on  an  axially  symmetric  body  and  a  flat  plate  is  found  to  be, 


Cf  /Cf\ 

T  \TJt.  P. 


i 

"si  f  (s) 


V?  (14) 


The  term  in  curly  brackets  in  the  above  equation  can  now  bo 
identified  as  a  combined  pressure  gradient  and  body  shape  parameter.  For 
brevity  this  is  given  the  symbol  Y  and  is  referred  to  simply  as  the 
pressure  gradient  parameter.  Consequently  we  have:- 

-m/m  +  1 

f  (s).  ds  ~]  _ VI  (15) 


The  values  of  the  constants  a  ,  p  ,  m,  and  B  are  as  follows i- 


Laminar  Flow 

Turbulent  Flow 

'  a  i  3-85 

!  p|  0.3 

5  m  !  1 

j  b|  5-18  §S  +  3 

3.269  * 

0.825 

0.25 

1.875  Tw/To+2. 25 

A  computer  programme  was  written  to  evaluate  ^  for  hemispheres 
and  hemispherical ly  blunted  cones.  In  all  cases  the  modified  Newtonian 
pressure  distribution  was  assumed  in  the  computation  of  ^  .  Also,  it 
was  necessary  to  take  a  mean  value  of  Tw/To  in  the  calculations  since 
the  theory  assumes  isothermal  wall  conditions. 
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APPENDIX  VII 

Estimation  of  Wail  Tomp§£atu£e,§ 


Providing  that  the  skin  can  be  regarded  as  thei'mally  thin,  Hill 
showed  in  Ref.  (32)  that  tho  wall  temperature  at  the  instant  tm  is 
related  to  that  at  the  previous  instant  tm_-|  f  by  the  equation 

(t„)b  -  ha  Trm  +  (hTr  -  hT„  -  2GT„  /  6  )„_,  - VII  (1) 

(2/6)  hm 

where  G  =  ps  Cs  T 

4 

and  r  =  radiation  rate  =  e  a  T 

w 

The  variation  of  specific  heat,  c,  with  temperature  was  included 
in  the  calculations. 

The  heat-transfer  coefficients  for  laminar  and  turbulent  flow 
were  determined,  respectively,  from  the  following  semi-empirical  relations 
(Ref.  33), 

h  =  0.0368  (Res)  1/5  (T-)  /Tw)°'46  (p  ^  Cp) 

'll 

(turbulent  flow)  _ VII  (2) 

h  =0.413  (Res)"^  (^/Tw)0-15  (puC) 

1  1  ! 

(laminar  flow)  _ VI  (3) 


where  a  Prandtl  n'umber  of  0.7  has  been  assumed. 


A  computer  programme  was  written  to  solve  the  above  set  of 
equations  with  the  iterations  repeated  until 


;w 


m 


-  T. 


1 


w;  m-1 


<  1 


+  0. 1°K 
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